Phase engineering of squeezed states and controlled entangled number states of 
Bose-Einstein condensates in multiple wells 
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We provide a scheme for the generation of entangled number states of Bose-Einstein condensates 
in multiple wells, and also provide a novel method for the creation of squeezed states without severe 
adiabatic constraints on barrier heights. The condensate ground state in a multiple well trap can be 
evolved, starting with specific initial phase difference between the neighboring wells, to a state with 
controllable entanglement. We propose a general formula for appropriate initial phase differences 
between the neighboring wells that is valid for any number of wells, even and odd. 
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Entanglement, a nonclassical correlation between two 
or more physical systems, lies at the heart of the pro- 
found difference between quantum mechanics and a lo- 
cal classical description of the world 0. Apart from 
their discussions in the philosophical and foundational 
aspects of quantum mechanics 2j, entangled states in 
recent years have become an essential resource for the 
emerging field of quantum information processing. En- 
tangled states have been created with photons [3j, four 
atoms 0, and most recently with many cold atoms in 
a Mott insulating state in an optical lattice |^. Cold 
atoms in optical lattices have been a vibrant research area 
with several new observations such as the superfluid to 
Mott insulator transition and number-squeezed states 
of Bose-Einstein condensate (BEG) 7]. Entangled and 
squeezed states hold promise in studies related to quan- 
tum measurement, the Heisenberg limited atom interfer- 
ometry and quantum computing and quantum commu- 
nication protocols 0- While the consequence of entan- 
glement for an Einstein-Podolsky-Rosen (EPR) pair is 
quantified in Bell's inequality 0, a more striking conflict 
between quantum mechanics and local realism is exhib- 
ited by three maximally entangled particles also known as 
the Greenberger-Horne-Zeilinger (GHZ) states jlJi]. GHZ 
state of N particles has the form 



entangled number states of a multiwell BEG of the form 
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where 1 and 2 are the basis states for a two state (spin 
^) system, and are written in the standard notation as 
|00...0) -I- In the occupation number basis, jl)"^^ 

and |2)'»^ respectively denote, \N,0) and \0,N). The 
two-state model has been generalized to more than two 
spin components in Ref. . The superposition of two 
macroscopicaiiy distinct states, rather than simply the 
internal degrees of freedom, each occupied by all N par- 
ticles, has been discussed by Schrodinger in the famous 
cat parable partial realization of such states has 

been obtained with Josephson junction loops (13) ■ 

In this paper, we discuss the generation of macroscopic 
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(|l)®^ + |2)®^ + ...-f|M)®^) (2) 



where 1, 2, 3, ..M label the macroscopicaiiy and spatially 
distinct wells, now denoting \0,0...,N, = N,...0). 

We show that states approximating the extreme entan- 
gled states of Eq. ^ may be generated in a controlled 
fashion by time evolution of appropriately phase im- 
printed ground states of a multiwell BEG with periodic 
boundary conditions for M=3 and 4. We show that the 
choice of initial barrier heights, which determine the ex- 
tent of ground state number squeezing, and the rate of 
barrier ramping can be used to control the extent of 
entanglement of the final states. We also show that 
fully fragmented states can be generated via natural time 
evolution from the ground state following certain initial 
phase offsets. The creation of such fragmented states 
through phase engineering and without the severe adi- 
abatic constraints on the rate of barrier height change 
provides an alternative to the current experimental ap- 
proaches 0, Finally, based on results obtained for 
two, three, and four well configurations, we conjecture a 
generalized formula, for M wells, for the phase offset be- 
tween neighborhing wells appropriate for the generation 
of number entangled states. 

We approximate the physics of a BEG in a multiwell 
potential by the Bose-Hubbard model Thus 

H = -J ^(a|ai+i + al+ifli) + ^ ein^ 

i i 

+ \uY,ndn^-l) (3) 



where 



is the number operator, J is the nearest 



neighbor tunneling term, U is the on-site energy, and 
is the energy offset of the ith lattice. To simplify a the- 
oretical study, we make a one parameter approximation 
of the tunneling and mean-field terms: U/J= l/e~"; 
and, for the symmetric wells explored here, Ci — Q. a 
is a dimensionless parameter that can be mapped onto 
the barrier height. This paramctrization allows a simple 
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FIG. 1: Fock state coefficients for 12 particles in three wells: 
(a) the ground state, (b) 10th, (c) 76th and (d) 91st, the 
highest state. The ground state has a Gaussian shape, while 
higher lying states are entangled number states, ni and n2 
are the Fock state indices and the vertical axis shows proba- 
bilities. Points beyond the cross-diagonal are unphysical. 



study of continuous change of barrier height through the 
variation of a single parameter a. For example, for a 
lattice made of red detuned laser with A ~ 985 nm and 
for ^•^Na, a barrier height IbE^j gives U — O.lbEn and 
J = O.OJEr Q where Er = is the recoil energy 
from absorption of a photon; these experimental param- 
eters then correspond to a = 2.14. 

In order to gain insight into the types of stationary 
states possible for the multiwell Bose-Hubbard model, 
we first analyze the quantum mechanical properties of 
the simplest multiwell potential, M = 3, assuming three 
symmetric wells in a circular array |l5l | . The state vector 
is a superposition of all the number states 
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Here ni, n2, and — N — ni — n2 are the number 
of particles in each of the three wells. Fig. ^ shows 

the Fock space probabilities, Cn[^n2 , for representative 

stationary states for = 12 and a = {U/J = 1). 
Our method of graphical representation is described in 
the figure caption. The ground state in Fig. ^a) is a 
broad Gaussian while the higher lying states, Figs.Q^b)- 
(d), are number entangled states of increasing extrem- 
ity corresponding to increasing numbers of particles si- 
multaneously in all three wells, the highest of which in 
panel (d) is an extreme superposition state of the form 
I A^, 0, 0) + |0, N, 0) + |0, 0, N). The number of non van- 
ishing Fock state coefficients determines sharpness, and 
thus (d) is sharper than (c). 

It is unlikely that such maximally entangled states can 




FIG. 2: Evolution to an entangled Fock space state: (a) bar- 
rier ramp showing the location of the following time evolved 
states: (b) initial state, (c) at 0.43 ms the Gaussian distri- 
bution broadens, (d) at 0.74 ms the distribution is 'splitting', 
(e) A three-peaked state is formed at 1.36 ms; a macroscopic 
superposition of definite number of particles simultaneously 
in all three wells. 



be generated via a sequence of single particle excitations. 
They may however, be dynamically generated via phase 
engineering from the approriate ground state. Writing 
phases on part of a condensate is experimentally feasible 
via interaction with a far off- resonance laser , and is 
assumed to be sudden with respect to the dynamics of 
the condensate . Mathematically, this corresponds to 
multiplying the coefficients in an expansion of the type of 
Eq.(4) by e™*^*, where jm, n2, ■■■rii, ...) is the correspond- 
ing Fock state, and 9i is the phase for particles in the ith 
well. Entangled state generation, obtained via integra- 
tion of the (linear) time-dependent Schrodinger equation, 
is shown in Fig.EJ following phase imprinting of an initial 
phase difference of ^ between the neighboring wells, and 
a simultaneous linear ramping of the barrier as a = i, as 
shown in Fig. [2Ia) (t here is dimensionless). Panels 3(b) 
shows the initial ground state; 3(c) at time 0.43 ms, the 
distribution broadens; 3(d) at 0.74 ms, in the process of 
splitting the state towards the three corners; and 3(e) at 
1.36 ms a sharp, although not extreme, entangled num- 
ber state with its signature of three major non vanishing 
expansion coefficients. The times are given for a *^Rb 
condensate, A = 840 nm, Ugc — 5.8 nm, J — 0.04i?fle~" 
and taking U — O.OAEn as approximately constant for 
calculational purposes. When an appropriately entan- 
gled state is reached the barrier is suddenly raised to halt 
further evolution in n-space. For the parameter values 
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FIG. 3: Entangled states evolved from ground states with 
different initial squeezing. Row (1) shows the states with 
a = 1.5 + t: (a) initial ground state and (b) final state. Row 
(2) is for a = 0.5 + t and (3) is for a = t. Column (b) gives the 
states at t=1.85 ms, t=1.36 ms, and t=0.99 ms respectively. 
The initial squeezing of the ground state thus determines the 
extremity of the resulting entangled states. 
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FIG. 4: Four well stationary and time evolved states: (a) 
ground state, (b) highest excited state, (c) phase engineered 
fragmented state following a 7r/2 relative phase shift of the 
ground state, (d) an entangled state evolved from the ground 
state following a tt relative phase shift. The three dimen- 
sions show the Fock state indices ni, n2 and na, probabilities 
are shown in the color intensity scale. For graphical clarity, 
only the points higher than 40% of the highest probability are 
shown, with the highest probabilities normalized to 1. 



used here, a simple time evolution without any change 
of barrier also produces an entangled state, however bar- 
rier ramping is used here to sharpen the resulting state. 
Control of the extremity of the states can be achieved by 
choice of the initial barrier height, controlling the initial 
squeezing of the ground state. This is demonstrated in 
Fig. 121 where different initial squeezing have been used 
for rows (1), (2) and (3). The columns show: (a) the 
ground state, and (b) the final state at the end of the 
barrier ramping. It is important to be able to tune to 
less extreme entangled states, as such states are more ro- 
bust to loss and decoherence |3 ■ Phase imprinting with 
a phase difference of ^ produces an equivalent state, 
with different phase space dynamics. 

The physics of the creation of these entangled states 
can be understood in terms of the underlying classical 
phase space dynamics. Entangled state generation in a 
double well has been thor oug hly analyzed in a semiclas- 
sical phase space picture [l7|. In the semiclassical limit 
valid for large N, the operators di can be approximated 
by the c-numbers y^e*^' , where rii and 6i are the num- 
ber and phase of particles in the ith well. The double 
well dynamics is then described by the Hamiltonian of 
a nonrigid physical pendulum [l^ with the number and 
phase differences (An,A0) between the wells as conju- 
gate variables. This system has two fixed points - (0,0) 
and (0,7r). The (0,0) is a stable equilibrium, while the 
(0,7r) is stable in the 7r-state regime {UN /J < 1) and 



unstable otherwise {UN/ J > 1). Exploitation of the bi- 
furcation characteristics of the unstable equilibrium gen- 
erates entanglement. Taking the initial parameters such 
that there is an unstable hyperbolic fixed point, phase 
imprinting moves the ground state to the unstable point, 
and the wave packet splits in the subsequent time evo- 
lution. Control of the barrier height can then be used 
in three different ways to control the motion of the wave 
packet and thereby the nature of the desired entangle- 
ment |l7j . First, a simultaneous ramping of the barrier 
with the natural dynamics at the unstable fixed point has 
been empirically found to be useful in directing the de- 
sired evolution of the wavepacket. Second, initial barrier 
height, that is the initial squeezing, helps shape the initial 
wave packet stretching it into different regions of acces- 
sible phase space; and, thirdly the initial barrier height 
sets the (negative) curvature of the potential at the hy- 
perbolic fixed point, controlling the rate of splitting of the 
wave packet. All of these effects can be visualized for the 
double well system in the appropriate phase space [l^ . 
Similar to the simple double well, the triple welL,_M = 3, 
can be thought of as two coupled pendulums [l9|] with 
complicated dynamics. In the {An,A9) representation, 
the unstable fixed points are {Ofi,^,^), k = 1,2. All 
the features of the double well entanglement generation 
apply to the three well case, and thus, many of the in- 
sights from the two and three well dynamics can be ex- 
tended to arbitrary number of wells in a circular array. 
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Next, we explore the dynamics of the ^ and tt phase 
shifting for the symmetric four weU case. First we look at 
the ground and the highest excited state in Figs.^a) and 
(b) for N = 16,U = 0.25, J = e"", a = 0.175, and as- 
suming the case of ®^Rb in the previous example. As ex- 
pected, the ground state is an approximate Gaussian and 
the highest state is an extreme entangled number state. 
The fixed point dynamics for the ^ configuration is such 
that, for a constant barrier height, the state evolves into 
a number-squeezed state during its evolution towards an 
entangled state. Fig.^fc) shows a fragmented state (at 14 
ms), with essentially exactly 4 particles in each well, ob- 
tained by this phase engineering scheme. Due to the loca- 
tion of the fixed points in the phase space, the entangled 
states generated by the ^ imprinting are not as extreme 
or sharp as those generated by the tt method; an example 
of the latter is shown in Fig. 2Jd) (at 1.5 ms). We also 
found that a ^ phase difference between the wells pro- 
duced entangled states that only differ by a symmetry 
from the ^ phase imprinted states. Fragmented states 
have been observed in a 12 well optical lattice with adia- 
batic raising of the barrier ,7] . We have shown here that 
fragmented states can also be created in a natural and 
efficient fashion without adiabatic constraints. In com- 
paring our results to the results of Ref . '2^ , we show that 
the TT configuration in an even number of wells that they 
identified is a special case of many phase imprint dynam- 
ics that generate interesting correlated states in multiple 
wells. Their changes in system parameters is to drive the 
system from stability to a regime of instability. On the 
other hand, we take our system to be in the unstable 
regime and demonstrate the controllability of entangled 
states with barrier manipulation; potentially useful for 
experimental detection. 

For the two, three, and four wells, we find M — 1 dis- 
tinct phase differences between the neighboring wells for 
the multiwell fixed points Ell | . These are given by a 
general formula where j = 1,2, ....M — 1, with M be- 
ing the number of wells, which gives a tt phase difference 
for the M — 2 double well, a ^ and ^ phase difference 
for the M = 3 triple well, and a ^, tt, and ^ phase 
difference for the M = 4 quadruple well configuration - 
we demonstrated the dynamics generated by all of these 
phase difference imprints. Note that the total change in 
phase in the circular loop is a multiple of 27r, a vortex 
like condition. We thus propose a general formula for M 
wells. 



A9 = 
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for the constant phase offset between neighboring wells 
leading to the dynamical generation of entangled states. 
Here j — 1,2,.., A/ — 1, and Eq. being valid for 
any number of wells, even or odd, provides a substantial 
generalization of the phase offset mentioned in Ref. , 
which is valid only for the special cases of an even number 



of wells and for j = M/2 (tt phase offset). The multiplic- 
ity of Eq. (O is prominent for large number of wells, e.g. 
for 12 wells, there are 11 phase offset possibilities. As in 
the three and four well case considered, symmetries may 
prevent all the imprinting offsets of Eq.(5) from generat- 
ing independent dynamics. 

In conclusion, we have demonstrated phase engineering 
schemes for the generation of entangled number states 
and fragmented states of BECs in multiple wells. By 
controlling the initial barrier height and rate of ramping, 
the entanglement of the final state can be tuned. We pre- 
sented a novel series of formulae for the initial phase dif- 
ference between the neighboring wells that is valid for any 
number of wells, even or odd, each having distinct proper- 
ties. The creation, characterization, and applications of 
multidimensional/multipositional Schrodinger cat states 
of atoms remain largely unexplored experimentally, and 
the theoretical ramifications of such states, should they 
be easily produced, are just emerging. 
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